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Abstract 

A possibility to prove spin and CP-eigenvalue of a Standard Model 
(SM) Higgs boson is presented. We exploit angular correlations in the 
subsequent decay H — ► ZZ — ► 41 (muons or electrons) for Higgs masses 
above 200 GeV. We compare the angular distributions of the leptons origi- 
nating from the SM Higgs with those resulting from decays of hypothetical 
particles with differing quantum numbers. We restrict our analysis to the 
use of the Atlas-detector which is one of two multi-purpose detectors at the 
upcoming 14TeV proton-proton-collider (LHC) at CERN. By applying a 
fast simulation of the Atlas detector it can be shown that these correla- 
tions will be measured sufficiently well that consistency with the spin-CP 
hypothesis 0+ of the Standard Model can be verified and the 0- and 1=L 
can be ruled out with an integrated luminosity of 100 fb . 
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1 Introduction 



Although the standard SU(2)l X U(l)y electroweak gauge theory successfully 
explains all current electroweak data, the mechanism of spontaneous symmetry 
breaking has been tested only partially. Since in the Standard Model sponta- 
neous symmetry breaking is due to the Higgs Mechanism, the search for a Higgs 
particle will be one of the main tasks of future colliders. 

At present, LEP gives a lower limit of 114.4 GeV/c 2 for the Higgs boson 
mass. There is also an indirect upper limit from electroweak precision measure- 
ments of 219 GeV/c 2 at a 95% confidence level [2] which is valid in the minimal 
standard model. However, this limit is still preliminary and the quality of the 
SM fit, when including all EW measurements from both low and high energy 
experiments, is still an object of discussion amongst experts Furthermore, a 
heavier Higgs boson would be consistent with the electroweak precision measure- 
ments in models more general than the minimal standard model |3]. In this first 
analysis we will therefore also consider higher Higgs masses well above this limit, 
as can be produced at high energy hadron colliders such as the LHC (pp collisions 
at 14 TeV). A Standard Model Higgs boson lying below the WW threshold will 
mainly decay into a bb pair. In this case, there is an overwhelming direct QCD 
background which dominates the signal. Therefore, the Higgs boson is difficult 
to study in detail in this mass region, even though one can use rare decays as a 
signal. Rare decays considered in the literature include, for example, H — ► t + t~, 
H — > 77, H — > Zj, H — > ZZ* or H — > WW*. All of these signals are rather 
difficult to see, but can eventually be used to establish the existence of the Higgs 
boson While these decay modes can be used to discover the Higgs boson, a 
detailed study of its properties will be difficult. 

The situation is much better for a heavy Higgs boson (m# > 2mw). For 
such a Higgs boson the main decay products are vector boson pairs, W + W~ or 
ZZ. For the latter decay mode, a clear signal for the Higgs consists of a peak in 
the invariant mass spectrum of the produced vector bosons. The double leptonic 
decay of the Z boson, H — > ZZ — ► Ifl^l^l^ , leads to a particularly clean signal. 

In this case, the basic strategy for discovering a Higgs boson in a clean mode 
is to select events with 4 high leptons that can be combined to form two 
Z-bosons. Here, an exposure of 30 fb _1 is already sufficient. If one finds such a 
signal one might be tempted to assume this to be the Standard Model Higgs bo- 
son. However, given the fact that the Higgs sector is not fully prescribed, one has 
to allow for other possibilities. In strongly interacting models, for instance, low 
lying (pseudo-)vector resonances are possible Cj. Also, pseudoscalar particles 
are present in a variety of models . Therefore, the first priority after finding a 
signal is to establish the nature of the resonance, in particular its spin and CP- 
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eigenvalue. This can be done by studying angular distributions and correlations 
among the decay leptons. In the following, we will make this study. We will 
limit ourselves to (pseudo-) vector and (pseudo-) scalar particles. To demonstrate 
consistency with a Spin 0, CP even hypothesis, we will compare the angular 
distributions, to those produced by different particles, always assuming the pro- 
duction rate of a Standard Model Higgs boson. This is the right assumption to 
make because, in order to be recognized as a candidate for a Standard Model like 
Higgs, the detected signal must be a resonance with the appropriate width and 
branching ratios. Since the production mechanism - gluon-fusion rules out spin 1 
particles, due to Yang's theorem[8j - cannot be seen, the only way to prove that 
the spin and CP nature of the new particle is Standard Model like is to study 
the decay angles of the leptons. 

Theoretical studies of angular distributions have been performed in the lit- 
erature [9-15]. So far, such studies have been limited to theoretical discussions. 
However, it was shown in (TJiJ that acceptance and efficiencies of the detector 
can play a role since they can generate correlations, mimicking physical ones. 
Therefore, it is necessary to use a detector simulation in order to establish how 
well one can do in practice. 

The complete triple differential cross-sections for a Higgs-boson decaying into 
two onshell Z-bosons which subsequently decay into fermion pairs can be calcu- 
lated at tree level. The angular dependence of this cross section is given in the 
appendix together with the most important integrated angular distributions. For 
the definition of the angles see Figure 1. 

We study essentially two distributions. One is the distribution of the cosine 
of the polar angle, cos 9, of the decay leptons relative to the Z boson. Because 
the heavy Higgs decays mainly into longitudinally polarised vector bosons the 
cross section da/dcos8 should show a maximum around cos(6 l )=0. The other is 
the distribution of the angle <fi between the decay planes of the two Z bosons in 
the rest frame of the Higgs boson. This distribution depends on the details of 
the Higgs decay mechanism. Within the Standard Model, a behaviour roughly 
like 1 + (3 cos 20 is expected. This last distribution is flattened in the decay 
chain H — > ZZ — > 41, because of the small vector coupling of the leptons, in 
contrast to the decay of the Higgs Boson into Ws or decay of the Z into quarks. 
Also, cuts can significantly affect the correlations. Therefore one needs a pre- 
cise measurement of the momenta of the outgoing leptons. The Atlas-detector 
should be well-suited to measure these distributions, since the muon and electron 
reconstruction is very precise over a large solid angle. A detector Monte-Carlo 
is however needed in order to determine whether the angular distributions can 
be measured sufficiently well in order to determine the quantum numbers of the 
Higgs particle. 
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Figure 1: The decay plane angle <p is measured between the two planes defined 
by the leptons from the decay of the two Z bosons in the rest frame of the Higgs, 
using the charge of the leptons to fix the orientation of the planes. The dashed 
lines represent the direction of motion of the leptons in the rest frame of the Z 
Boson from which they originate. The angles Q\ and 82 are measured between 
the negatively charged leptons and the direction of motion of the corresponding 
Z in the Higgs boson rest frame. 0=0 correspond to p e + x p e - and p M + x p^- 
being parallel. 0=7r correspond to p e + x p e - and p^+ x p^~ being antiparallel. 



The present paper is organized as follows. In Chapter 2 the generator is 
described, in Chapter 3 detector simulation and reconstruction are given. In 
Chapter 4 we define quantities that can be used to characterize the different 
distributions. In Chapter 5 we present the results, concluding that the quantum 
numbers of the Higgs particle can indeed be determined. In the appendix we give 
formulae for the complete differential and integrated distributions for the decay 
of the resonance assuming arbitrary couplings computed in tree level and narrow 
width approximation. 

2 The Generators 

In order to distinguish between different spins J=0,1 and/or CP-eigenvalues jcp — 
—1,-1-1 one needs to study four different distributions: that resulting from the 
decay of the Standard Model Higgs boson, and the three distributions that would 
result from hypothetical particles with spin and CP-eigenvalue combinations (0, 
1), (1, 1), (1, -1). 

The feasibility of using angular correlations in the decay of the Z bosons in 
order to distinguish between these particles has been evaluated using two dif- 
ferent Monte-Carlo generators. One was written for the Standard Model Higgs 
(gg — > H — > ZZ — > 41) and the irreducible ZZ-background. 14j The latter includes 
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contributions from both gluons and quarks to the ZZ production ( gg — > ZZ — > 4/ 
and —* ZZ —* 41) , whereas all Higgs production mechanisms other than the 
gluon fusion are neglected. This generator keeps all polarisations of the Z-boson 
for the quark initiated as well as for the gluon initiated processes. This allows for 
an analysis of the angular distributions of the leptons. The gluonic production of 
Z-boson pairs is only about 30% of the total background. However, one should 
not ignore its contribution, since it has different angular distributions from the 
other backgrounds and its presence can affect measured correlations. The pro- 
gramme contains no K-factors, therefore our conclusions regarding the feasibility 
of the determination of Higgs quantum numbers are conservative. Indeed, K- 
factors are expected to be larger for gluon-induced processes (such as the signal) 
than for quark-induced processes (70% of the background). Since the narrow 
width approximation is used, the results are only valid for Higgs masses above 
the ZZ threshold. 

For the alternative particles, a new generator was written based on an article 
by C. A. Nelson and J. R. Dell' Aquilla ^3]- The programmes for the produc- 
tion of background, the Standard Model Higgs and all alternative particles use 
Cteq4M structure functions ^Hj and hdecay ^Jj for branching-ratio and width of 
the Higgs, and all use the narrow width approximation. The background as well 
as all cross sections for the four simulations are taken from the first generator. 
Thus, all cases show identical distributions of invariant mass of the Z-pairs and 
transverse momentum Pj< of Z-bosons and leptons and have the same width and 
cross sections. The only difference lies in the angular distributions of the leptons. 
For the alternative particles no special assumption concerning the coupling has to 
be made; only CP-invariance is assumed. It is worth mentioning that the angular 
correlations are completely independent of the production mechanism. 



3 Detector Simulation and Reconstruction 

The detector response is simulated using ATLFast [TH], a software-package for 
particle level simulation of the Atlas detector. It is used for fast event-simulation 
including the most crucial detector aspects. Starting from a list of particles in the 
event, it provides a list of reconstructed jets, isolated leptons and photons and 
expected missing transverse energy. It applies momentum- and energy- smearing 
to all reconstructed particles. The values of the detector-dependent parameters 
are chosen to match the expected performance that was evaluated mostly by full 
simulations using Geant3[19j. 

The event selection is modeled exactly after the event selection in the Atlas- 
Physics-TDR [20]. Four leptons (electrons or muons) are required in the pseu- 
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dorapidity range \rj\ = | lntan( B ^ am )| < 2.5 {9 Beam being the angle to the beam 
axis). Two of the leptons are required to have transverse momenta greater than 
20 GeV/c and the two other leptons must have transverse momenta greater than 
7 GeV/c each. A lepton identification efficiency of 90% per lepton was assumed. 
Two Z bosons are reconstructed by choosing lepton-pairs of matching flavour 
and opposite sign. If the flavours of all four leptons are equal, the combination is 
chosen, which minimizes the sum of the squared deviation of the invariant mass 
of the pairs with opposite sign from the Z mass (i.e. choose combination ab/cd 
that minimizes (m a +6--m^) 2 + (m c + a —-rciz) 2 )- The reconstructed invariant mass 
of the two reconstructed Z bosons has to lie within two times the width of the 
reconstructed mass peak of the Higgs resonance around the centre of the peak. 
For high Higgs masses (m# > 300 GeV/c 2 ) this is only little more than two times 
the decay width, while for smaller masses the experimental resolution dominates. 

Throughout this paper, we use the term signal for distributions where the 
background has been statistically subtracted. The only background considered 
is the Z pair production. Other possible backgrounds like top pair production 
or Zbb are negligible for masses of the Higgs boson above 200 GeV/c 2 . System- 
atic uncertainties due to the simulation of the background could be studied by 
comparing distributions from the sidebins of the Higgs signal with the results of 
the generator. A proper treatment of the background is very important, since 
the angular distributions of the background itself and correlations introduced by 
detector effects have a large impact on the shape of the distributions discussed. 
These effects are detailed below. 

For high invariant masses, the Z bosons from the background processes are 
mainly transversely polarised leading to a polar angle distribution of the form 
T^-g ~ 1 + cos 2 9. This distribution flattens the sin 2 9 distribution expected for 
the Higgs decay. Figure El shows the polar angle distributions of the signal (left) 
and the background (right). The dashed line shows the shape of the distribution 
expected when no cuts are applied and the detector response is not taken into 
account. It has just been scaled by the overall acceptance of the cuts, so that 
the shape can be compared. The expected distribution with all cuts and smear- 
ing applied is drawn as a solid line. Figures El and 01 are produced assuming a 
Higgs mass of 200GeV/c 2 and Z decaying to muons only. For the decay ZZ — > 
e + e~e + e~ or ZZ — > /i + /i~e + e~ the graphs look similar. The smearing effects are 
largely independent of the Higgs mass. 

The effect of the detector acceptance and isolation cuts on the decay plane 
angle distribution is shown in Figure El Again, the distribution of the signal 
is shown left and the background right. The dashed histogrammes are scaled to 
have the same integral as the solid histogrammes, and zero is suppresssed in order 
to facilitate the comparison of the shape. The definition of the line styles are the 
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Polar angle distributions of signal and background 
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Figure 2: Distribution of the polar angle cos(#) for the background only (right) 
and the signal (left). The Higgs mass is 200 GeV/c 2 . 
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Decay plane angle distributions of signal and background 
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Figure 3: Distribution of the decay plane angle for the background only (right) 
and the signal (left). The Higgs mass is 200 GeV/c 2 . The distribution without 
cuts is scaled by the expected acceptance, so that the shape can be directly 
compared. 



same as above. For the decay plane angle, the background shows an almost flat 
distribution before applying selection cuts. But a minor correlation is introduced 
by detector effects. This has been simulated and taken into account for the anal- 
ysis of decay plane angle distributions. In conclusion, the isolation cuts lead to a 
small distortion of the angular distributions as discussed in [TJ] , but these effects 
are almost negligible for the Atlas detector. The cut on \rj\ enhances the decay 
plane correlation a little, but the smearing and the Pt requirements reduce this 
effect. Altogether, there is a small enhancement of the correlation of almost the 
same amount for all four particles. 

A further cut on the transverse momentum of the Z bosons P T nax {Zi, Z 2 ) > 
ran 1 3 is known to additionally reduce the background, but it also affects the 
correlation. Since an optimisation of the signal-to-background ratio is not crucial 
to this analysis, this cut has not been applied, rendering the analysis less depen- 
dent on the details of the production mechanism like initial Pt of the Higgs boson. 
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Table 1: Distribution of the polar angle 9. Pi are the Legendre Polynomials. See 
the text for definitions. 



4 Parametrisation of Decay Angle Distributions 

The differential cross sections for the different models can be computed directly 
or can be derived from the formulae given in [T^j. The explicit distributions are 
given in the appendix. From the article [12, we quote the simple distributions 
of the alternative particles. Table d shows the distribution of the polar angle 9. 
9\ and 9 2 are the polar angles of the leptons originating from the Z Bosons Z\ 
and Z 2 respectively. In Table EJ the distribution of the decay plane angle <fi is 
shown where the polar angle 9 is integrated over different ranges. Fll gives the 
distribution for < lj2 < tt/2, F22 for tt/2 < 9 1}2 < tt, F12 for < 9 1 < tt/2 
and tt/2 < 9 2 < tt, and F21 for tt/2 < 9 X < tt and < 9 2 < tt/2. K, U, T and 
W are the parameters that characterise the decay density matrix. For the decay 
modes used in this analysis, they amount to the following values: 1Z = U = -1/2, 
T = W = — jT^r- r is the ratio of the axial vector to vector coupling which for 
the muons amounts to r = (1 - 4 sin 2 We used sin 2 ^ = 0.23. 



The plane-correlation can be parametrised as 

P(0) = 1 + a ■ cos(0) + (3 ■ cos(20) (1) 

In all four cases discussed here, there is no sin(0) or sin(2</>) contribution. For 
the Standard Model Higgs, a and (3 depend on the Higgs mass while they are 
constant over the whole mass range in the other cases. 



The polar angle distribution can be described by 

G{9) = T- (1 + cos 2 (fl)) + L-sin 2 (#) (2) 

reflecting the longitudinal or transverse polarisations of the Z boson. We define 
the ratio 

R--=^- (3) 

L + T y ' 

of transversal and longitudinal polarisation. 

The dependence of the parameters a, (3 and R on the Higgs mass is shown in 
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Figure 4: The variation of the three parameters a, (3 and R (top to bottom) with 
the Higgs mass. 
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Table 2: Distribution of the decay plane angle 0. Fll gives the distribution for 
< #1,2 < 7r/2, F22 for tt/2 < #i >2 < vr, F12 for < 6\ < vr/2 and tt/2 < # 2 < vr 
F21 for tt/2 < ^ < vr and < 6 2 < tt/2. TZ = U = -1/2 , T = W = -jf^, r = 

(l-4sin 2 ^) -1 . 

Figure |U The pseudoscalar shows the largest deviation from the SM Higgs. It 
would have j3 = —0.25 and R — — 1 whereas the scalar always has (3 > and 
R > 0. The vector and the axialvector can be excluded through the parameter 
R for most of the mass range, but for Higgs masses around 200 GeV/c 2 the main 
difference lies in the value of j3 which is zero for J = 1 and 7c*p = ±1 and about 
0.1 for the scalar. The value for a can only discriminate between the scalar and 
the axialvector but the difference is very small. 



5 Background Estimation 

The subtraction of the angular distributions of the background is necessary to 
obtain and analize the angular distributions of the signal alone. This bears a 
risk of introducing systematic errors. Thus, the background distributions as pro- 
duced by Monte-Carlo-Generators have to be checked against the data. In this 
chapter we will estimate the effects and possible systematic errors introduced by 
the subtraction. 

First, the absolute number of background events has to be estimated. This 
can be done by comparing the sidebands of the signal to a simulated distribution 
of the background only. This procedure is illustrated in Figure El In order to ob- 
tain the number of expected events the number of simulated events in the signal 
region N*f^ ai is scaled by the number of events in the sidebands N^; a divided 
by the number of simulated events in the sidebands N^£. The error from this 

calculation is cr N = ^JN^T • -j^wu 1 - In the case of a 250 GeV/c 2 Higgs boson 
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4 Lepton Invariant Mass Distribution 




Figure 5: The invariant mass distribution of a 250 GeV/c 2 Higgs boson and the 
ZZ Background. The vertically hatched region is the signal region used in the 
analysis. The diagonally hatched regions are the sidebands used to determine the 
expected number of background events (hatched horizontally) inside the signal 
region. The dotted line indicates the shape of the background in the transition 
region between the sidebands and signal which is not used at all. 
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Measurement of Parameter R j 
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Figure 6: The parameter R (defined in (3)) as obtained by a fit to different mass 
regions of the background only (solid line) and by a fit to the same mass regions 
to signal plus background distributions (points with errorbars). The horizontal 
errorbars indicate the regions from which the distributions where taken. 
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the estimated number of background events is N=130 with a systematic error of 
CNsyst = 4.1 which is well below the statistical error of ciNstat = 11-4. 

Checking of the shape of the background distribution can be done by using 
bins below and above the signal region, too. This is demonstrated in Figure 
El It shows the parameter R derived from a fit described in chapter 0] to the 
background distribution only (black line) and the background plus signal (points 
with errorbars). For most of the fitted values a bin width of 20 GeV/c 2 was 
used, except for the last bin where 100 GeV/c 2 were used to compensate for the 
fact that there are less events for higher invariant masses. From the expected 
errors one finds that the parameter R for the background can be estimated with 
a precision of about ctr = 0.08. This might not seem too good, but the effect of 
using a slightly wrong background distribution is not so large. To demonstrate 
this, a fit to the angular distribution of the angle 9 was performed, where a 
wrong background distribution was subtracted from the signal-plus-background 
distribution as obtained from the generator. The parameter R su b of the subtracted 
distribution was changed to values higher and lower than the value of Kmc of 
the generated distributions. In Table El the difference AR = Rmc - Rs«ft and 
the value of R S ignai obtained from the fit to the signal distribution produced by 
subtracting the wrong background distribution are shown. 



AR 
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-0.1 
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Rsignal 
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0.770 


0.782 


0.796 



Table 3: The measured Parameter R for five different distributions that have 
been used to subtract the expected background distribution. AR is the difference 
between the value of R from the background as produced by the Monte Carlo 
and the value of R of the subtracted distribution. 

The shift of the parameter is thus expected to be less than ±0.01. Again, 
this is very small compared to the statistical error of AR stat = 0.053. This error 
is not considered in the rest of the analysis. Furthermore, the effects will be 
even smaller when considering K-factors. Any K greater than 1 will give better 
conditions to check the background distributions. And, since the K-factor of the 
gluonic Higgs production is higher than the K-factor of the main ZZ production 
process by quark antiquark pairs, the signal to background ratio will be even 
higher than predicted here. 

6 Results 

In Chapter 4, the exact results for the signal were given. However, in practice 
one needs a procedure to separate signal from background, which will lead to 
uncertainties in the distributions. The expected errors have been calculated by 
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Polarisation of the Z Bosons from Higgs decay (100 fb" 1 ) 



T± 1 



ii 

DC 

1 0.5 
E 

CO 
a- 
CO 
Q. 



-0.5 



if 



. SM-Higgs 
\j Spin 1, CP ±1 



Spin 0, CP -1 



_l I I L_ 



200 



250 



300 



m H [GeV] 



Parameter a and p 100 fb" 1 m H = 200 GeV ( 196 < l\<| nv < 204 ) | 



0.2 
0.1 

-0.1 
-0.2 
-0.3 













SM - Higgs 


Spin 1, CP -1 




Spin 1, CP 1 


i , i 


Spin 0, CP -1 

i , i , i 


-0.1 -0.05 







0.05 0.1 



a 



Figure 7: The parameter R for different Higgs massses (top) and a and (3 (bot- 
tom) for m H = 200 GeV/c 2 using 100 ftr 1 . The error scales with the integrated 
luminosity as expected. 



15 



Parameter a and p 100 fb 1 m H = 200 GeV ( 196 < l\^ nv < 204 )| 



CO. 



0.2 
0.1 
-0 
-0.1 
-0.2 
-0.3 
-0.4 



Spin 1, CP 1 



SM-Higgs 



Spin 1, CP -1 



Spin 0, CP -1 



_L 



-0.2 



-0.1 



0.1 



0.2 



a 



Parameter a and p 100 fb 1 m H = 200 GeV ( 196 < l\/| nv < 204 ) | 

co. 0.3 r 



0.2 
0.1 
-0 
-0.1 
-0.2 
-0.3 
-0.4 



SM - Higgs 



Spin 1, CP -1 



Spin 1, CP 1 



Spin 0, CP -1 



-0.3 



-0.2 



-0.1 



0.1 



0.2 



0.3 



a 



Figure 8: The parameter a depends on the signs of the cos(#) of the two Z bosons. 
The events where the signs are equal are used for the upper plot, those where the 
signs are different are used for the lower plot. 
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generating a large number of events and scaling the distributions to the expected 
number of events, since the expected values of the parameters follow a Gaus- 
sian distribution. The background was statistically subtracted after applying the 
same cuts to it as were applied to the signal. The error reflects the statistical 
error from the number of the signal events, the statistical error from the num- 
ber of background events subtracted and the error made by the estimation of the 
number of background events as described in Chapter 5. No error from a possibly 
different angular distribution of background events has been taken into acount, 
but we have shown that the effect is small. Then the parametrisations for and 
6 as described above were fitted to the distributions. Signal and background are 
summed over muons and electrons. 

Figure (top) shows the expected values and errors for the parameter R, 
using an integrated luminosity of 100 fb _1 . It is clearly visible that for masses 
above 250 GeV/c 2 the measurement of this parameter allows the various hy- 
potheses considered here for the spin and CP-state of the "Higgs Boson" to be 
unambiguously separated. 

For a Higgs mass of 200 GeV/c 2 only the pseudoscalar is excluded. Figure 
[3 (bottom) shows the expected values and errors for a and (3 for a 200 GeV/c 2 
Higgs and an integrated luminosity of 100 fb _1 . 

The parameter a can be used to distinguish between a spin 1 and the SM 
Higgs particle, but its use is statistically limited. The same applies to the pa- 
rameter (3. Measuring (3, which is zero for spin 1 and > in the SM case, can 
contribute only very little to the spin measurement even if m# is in the range 
where (3, in the SM case, is close to its maximum value. Nevertheless, (3 can be 
useful to rule out a CP odd spin particle. 

The values of a get more widely separated when the correlation between the 
sign of cos(6 l ) for the two Z Bosons and (ft is exploited. In Figure |H1 we plot 
the parameters separately for sign(cos9i) = signicosOz) (Fll + F22 in Tabled 
and sign(cos6i) = —sign{cos92) (F12 + F21 in Table HJ). As can be seen, the 
difference in a becomes bigger for J = 1 and 7cp = +1. For higher masses a 
and (3 of the SM Higgs approach 0; thus only a can be used to measure the spin. 
But the measurement of R compensates this. 

Figure El shows the significance, i. e. the difference of the expected values 
divided by the expected error of the SM Higgs. We add up the significance for 
a and j3 exploiting the cos(6 l ) - <fi correlation and plot the significance from the 
polar angle measurement separately. For higher Higgs masses the decay plane 
angle correlation contributes almost nothing, but the polarisation leads to a good 
measurement of the parameters spin and CP-eigenvalue. For full luminosity (300 
fb _1 ) the significance can simply be multiplied by V3. This is especially interest- 
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ing for a Higgs mass of 200 GeV/c 2 . The Spin 1, CP even hypothesis can then 
be ruled out with a significance of 6.4a, while for the Spin 1, CP odd case the 
significance is still only 3.9a. 

In conclusion, for Higgs masses larger than about 230 GeV/c 2 a Spin 1 hypoth- 
esis can be clearly ruled out already with 100 fb -1 . For m H around 200GeV/c 2 
the distinction is less clear, and one will need the full integrated luminosity of 
the LHC. A spin-CP hypothesis of 0- can be ruled out with less than 100 fb _ l 
for the whole mass range above and around 200 GeV/c 2 . 
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Figure 9: The overall significance for the exclusion of the non standard spin 
and CP-eigenvalue. The significance from the polar angle measurement and the 
decay-plane-correlation are plotted separately. 
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A Formulae for differential angular distributions 



The most general coupling of a (pseudo) scalar Higgs boson to two on-shell Z- 
bosons is of the following form: 

£>scalar = ^^fj,u + YA^A^/M^ + iPe fiUpzqz /Ml (4) 

Here the momentum of the first Z-boson is p z , that of the second Z-boson is 
q u z . The momentum of the Higgs boson is k and e^ vpcT is the total antisymmetric 
tensor with €1234 = %■ Within the Standard Model one has X = 1, Y = P = 0. 
For a pure pseudoscalar particle one has P^0,X = Y = 0. If both P and one 
of the other interactions are present, one cannot assign a definite parity to the 
Higgs boson. 

The same formula for a (pseudo) vector with momentum k p reads: 

^vector = ~X-(5 pll p z + S pu q z ) + P(ie /Ji ^p pz — it^upqz) (5) 

It is to be noted that the coupling to the vector field actually contains only two 
parameters and is therefore simpler than to the scalar. 

In the following we give the angular dependence of the triple differential cross 
section for the case of a scalar or vector Higgs decaying into two on-shell Z bosons 
which subsequently decay into two lepton pairs. The meaning of the angles 6\, 82 
and is explained in Figure ^ P is the absolute value of the momentum of the 
Z boson, p 2 = (^Mh) 2 — M|. In the following we use the definitions x = and 
y = -j^. c v and c a are the vector and axial vector couplings: c v = £3 — 2q sm(9w) , 
c a = £3, where £3 is the weak isospin, q the charge of the fermion and 0w the 
Weinberg angle. For our case, the values of c v and c a are c v = -0.0379 and c a = 
-0.5014. 
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A.l The general case 

Scalar Higgs 

da 

d(j)d cos 6 \d cos 2 

- 8XYc 2 a c 2 v x 2 (x 2 - 4) cos sin 9 X sin 2 

- XY(c 2 + c 2 a ) 2 x 2 (x 2 -4) (2 cos sin X sin0 2 cos0i cos0 2 + - 2) sin 2 X sin 2 2 ) 
+ lQXFc 2 a c 2 v xy (x 2 - 2) sin sin X sin 2 

+ 4XP(c 2 + c 2 ) 2 xy sin 0(2 cos sin 2 X sin 2 2 + {x 2 - 2) sin 9i sin 2 cos 6 1 cos 2 ) 

+ 16X 2 c 2 c 2 x 2 (2 cos 0i cos 2 + (x 2 - 2) cos sin X sin 2 ) 

+ X 2 (c 2 + c 2 ) 2 a; 2 {4(l + cos 2 Q x cos 2 2 + cos 2 sin 2 X sin 2 2 

+ (x 2 — 2) cos sin 0i sin0 2 cos0i cos0 2 ) + x 2 (x 2 — 4) sin 2 9\ sin 2 2 } 

- SPY c 2 a c 2 v xy(x 2 - 4) sin sin X sin 2 

- 2PY(c 2 + c 2 ) 2 xy(x 2 — 4) sin sin 0i sin 2 cos X cos 2 
+ l/4Y 2 (c 2 + c 2 a ) 2 x 2 (x 2 - 4) 2 sin 2 G 1 sin 2 2 

+ 8P 2 c 2 c 2 (a; 2 - 4) cos 9 1 cos 2 

+ P 2 (c 2 + c 2 a ) 2 (x 2 - 4)(1 + cos 2 0i cos 2 2 - cos 2 sin 2 9 1 sin 2 2 ) 

Vector Higgs 

da 

d(j)d cos 9 \d cos 9 2 

— 16XPc 2 c 2 x|/ sin sin 9\ sin 2 

+ 4XP(c 2 + c 2 ) 2 xy sin sin 9\ sin 2 cos 0i cos 2 
+ 4X 2 c 2 c 2 a; 2 cos sin 0i sin 2 

+ X 2 (c 2 + c 2 a ) 2 x 2 (l - cos 2 0i cos 2 2 - cos sin 0i sin0 2 cos0i cos0 2 ) 

- 4P 2 c 2 c 2 [x 2 - 4) cos sin 9 1 sin 2 

+ P 2 (c 2 + c 2 a ) 2 (x 2 - 4) (1 - cos 2 0i cos 2 2 + cos sin 9 X sin 2 cos 9 X cos 2 ) 
A. 2 The special cases 

In this appendix we list the triple differential cross section for pure Higgs spin and 
CP states. In addition, we also give the differential cross sections, where some 
of the angular variables have been integrated over. Fll, F12, F21, F22 refer to 
the different quadrants as defined in Chapter 4. The spin 0, CP even part only 
contains the pure SM contribution. 
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Spin 0, CP even 

— — — ~ + 16c 2 c 2 (2 cos 9i cos 9 2 + (x 2 - 2) cos0sin#i sin# 2 ) 

dcpd cos Bid cos 9 2 

+ ( c l + c a) 2 { 4 (! + cos2 81 cos2 #2 + cos 2 sin 2 9 l sin 2 6> 2 

+ (x 2 - 2) cos sin #1 sin# 2 cos#i cos# 2 ) + x 2 (x 2 - 4) sin 2 9 1 sin 2 9 2 } 

do 22 

— — ~ + 32c a c v cos 61 cos 9 2 

d cos &i d cos fl 2 

+ (c 2 + c 2 ) 2 {4(l + cos 2 9 1 cos 2 2 ) + (x A - Ax 2 + 2) sin 2 6 l sin 2 # 2 } 

Fll = F22: ^ ~c 2 c 2 (8 + 7r 2 (a; 2 -2)cos0) 
d<p 

+ 4/9 (c 2 + c 2 a ) 2 (x 4 - Ax 2 + 10 + (x 2 - 2) cos + 4 cos 2 0) 

F12 = F21: ^ ~ - c 2 c 2 (8 - tt 2 (x 2 - 2)cos0) 
a0 

+ 4/9(c 2 + c 2 a f{x A - Ax 2 + 10 - (x 2 - 2) cos + 4 cos 2 0) 



Spin 0, CP odd 

da 

d(j)d cos^icicos 6* 2 



da 



cicos 9±d cos 9 2 



+ 8c 2 c 2 cos #1 cos # 2 

+ (c 2 + c 2 ) 2 (l + cos 2 ^ cos 2 6» 2 - cos 2 sin 2 9 l sin 2 9 2 ) 

+ 16c 2 c 2 cos ^ cos 9 2 + (c 2 + c 2 ) 2 (l + cos 2 6>i)(l + cos 2 6> 2 ) 



Fll = F22: g ~ c 2 c 2 + l/9(c 2 + c 2 ) 2 (5 - 2 cos 2 0) 



F12 = F21: ^ ~ - c 2 c 2 + l/9(c 2 + c 2 ) 2 (5 - 2 cos 2 0) 



dcj) 
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Spin 1, CP even 

da 

d(j)d cos^icicos 9 2 



+ 4c 2 c 2 cos sin 6>i sin 9 2 

+ (c 2 + c 2 ) 2 (l — cos 2 #i cos 2 9 2 — cos sin 9± sin # 2 cos 9± cos 2 ) 



(icos^irfcos 9 2 



~ 1 - cos #i cos 9 2 



Fll = F22: ^ ~ + c 2 c 2 7r 2 cos0 + l/9(c 2 + c 2 ) 2 (32 - 4cos0) 
d(p 



F12 = F21: ^ ~ +c 2 c 2 tt 2 cos0+ l/9(c 2 + c 2 ) 2 (32 + 4cos0) 
d(p 

Spin 1, CP odd 

da 

r 

d(f)d cos 9 \d cos # 2 



— 4c 2 c 2 cos sin 9\ sin # 2 
+ (c 2 + c 2 ) 2 (l — cos 2 9\ cos 2 9 2 + cos sin 9\ sin # 2 cos 9\ cos # 2 ) 



rfcos^irfcos^ 



~ 1 - cos 9 1 cos 9 2 



Fll = F22: ^ c 2 c 2 tt 2 cos0+ l/9(c 2 + c 2 ) 2 )(32 + 4cos0) 

d<p 



F12 = F21: ^ c 2 c 2 tt 2 cos + l/9(c 2 + c 2 ) 2 ) (32 - 4 cos 0) 



d0 
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